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The physical principles are studied for the optimal design of a quantum multiwell heterostructure 
working as an electrostatic energy storage device. We performed the search for an optimal multiwell 
trapping potential for electrons that results in the maximum static palarizability of the system. The 
response of the heterostructure is modeled quantum mechanically using nonlocal linear response 
theory. Three main design strategies are identified, which lead to the maximization of the stored 
energy. We found that the efficiency of each strategy crucially depends on the temperature and 
the broadening of electron levels. The energy density for optimized heterostructures can exceed the 
nonoptimized value by a factor more than 400. These findings provide a basis for the development 
of new nanoscale capacitors with high energy density storage capabilities. 



Attaining high density electromagnetic energy storage 
is a vital technological and scientific problem, since elec- 
tricity is the most universal and scalable form of en- 
ergy. In addition, the continued progress in chip de- 
sign also produces demands for nanoscale capacitors with 
the enhanced energy density. Recently, different types of 
nanostructures, such as multi- walled doped carbon nan- 
otubes [H or metal-dielectric interfaces @, 0] have at- 
tracted much attention as potential building blocks for 
high energy density storage devices. 

There is no clear understanding of the physical prin- 
ciples for optimal design on scales where quantum me- 
chanical effects (e. g., like discreteness of energy levels, 
tunneling and nonlocality of dielectric response) start to 
play a significant role. A recent study [4( predicted the 
surprising effect that there is an optimal size for small 
metallic dimers resulting in maximal electromagnetic en- 
ergy density between the particles. This behavior is con- 
trary to the classical theory produces a monotonous in- 
crease of the energy density with decreasing size of a 
dimer. Similar simulations, but using more sophisticated 
ab initio techniques [f| confirmed the correctness of qual- 
itative predictions made in 0]. This confirmation pro- 
vides a foundation for performing a systematic study of 
the principles of optimal design of nanocapacitors using 
the formalism developed Q , thereby avoiding much more 
computationally expensive ab initio simulations. 

The present study focuses on the static limit of the 
dielectric response of heterostructures. A dynamic field 
may result in higher energy densities due to the reso- 
nance phenomena, however, this benefit is offset by much 
higher losses in the dielectric, which is not acceptable for 
energy storage applications. In the static regime, there 
are still energy losses due to electron-phonon coupling, 
etc. Typical values for the dissipation factor (reciprocal 
to the quality Q factor) is of the order 10 -3 or less. In 
the present study we assume zero dissipation factor in 
the static limit. 



This Letter considers a quantum multi-well het- 
erostructure working as an electrostatic energy storage 
device. We assume that the heterostructure is placed 
between two leads, without direct contact to the leads, 
permitting the neglect of tunneling effects. If a constant 
bias is applied to the leads, the electron density in the 
multi-well heterostructure becomes perturbed resulting 
in a finite polarization of the nanostructure. The en- 
ergy stored in the heterostructure depends on the di- 
electric function of the heterostructure, which is deter- 
mined by the eigenfunctions and eigenenergies of the con- 
fined electrons. The goal of optimal design is to find 
a shape for the heterostructure, which maximizes the 
stored electrostatic energy. In the present simulations 
a multi-well heterostructure is modeled by the Hamilto- 
nian H = — V 2 + Vt rap (r), where Vt rap (r) is the electron 
effective trapping potential. 

Since inhomogencous systems with a few electrons can 
have an electromagnetic response very different from the 
bulk, in our simulations we use the non-local density- 
density response function within the coordinate-space 
representation %(r, r', uj) : 



X(r,r',u) = ^2 



f(Ej) - f(Ej) 
Ei — Ej + hui + «7 



^(r)^(r')V*(r')^(r), 



(1) 



where f(Ei) is the Fermi filling factor and the small con- 
stant 7 describes the level broadening, and Ei and i\)i 
are the eigenenergies and eigenfunctions of the Hamilto- 
nian H, and u) is the frequency of the external bias field. 
The electron eigenenergies Ei and eigenfunctions tpi are 
obtained numerically using spacial discretization of the 
Schrodinger equation. Eq.JT]) permits calculation of the 
induced electric field Ei ra d(r, uj) in the system within lin- 
ear response theory jy]. 

The effective trapping potential Vt r a P (r) is represented 
as a sum of N p model potentials V mo d(r) located at 
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FIG. 1: The eigenspectrum for a non-optimized potential (a) 
and for the optimized one (b). The occupied states colored 
in blue and unoccupied states are shown in magenta. Note, 
that the energy gap AEn,n+i in (b) is much smaller than in 
(a), (c) The re-scaled trapping potential (black dashed line) 
and the corresponding ground state electron density (red solid 
line) for non-optimized and optimized (d) configurations, (e) 
The induced field intensity in the optimized trapping potential 
(red dashed line) is enhanced by several orders of magnitude, 
compared to the non-optimized potential (black solid line), 
(f) Frequency dependence of the induced field energy in the 
nanostructure. Note, T = 0, 7 = 0.01-En. 



the positions i = l,..,N p , which may be controlled 
during the material fabrication procedure: Vt rap (r) = 

Yn=i Vw(r - Ti). For simplicity we assume V mod {v) to 
be the Coulomb potential of a charge Q with a cutoff 
a: V mod {r) = Q/\r\, for |r| > a, V mod (v) = Q/a, for 
r < a. The algorithm searches for optimal parameters 
Yi to maximize the stored electrostatic energy in the sys- 
tem W = J v \E in d{r,uj — * 0)\ 2 dV 7 where the integration 
is performed over the volume of the nanostructure. The 
search for the optimal potential Vt rap (r) is performed us- 
ing Brent's "principal axis" optimization algorithm [7J 
for a scalar function of several variables. 

To simplify our simulations we consider a quasi-one 
dimensional system of length L. Then Eq = 



the unit of energy determined by the system, where m e 
is the electron mass. The optimization of the trapping 
potential is performed at the static limit of the external 
field hid = 0, with N = 15 electrons trapped in the het- 
erostructure having N p = 5 wells. The cut-off parameter 



is chosen as a — 0.05L, and Q — 3|e to ensure the total 
neutrality of the heterostructure. Note, that in all the 
simulations we assume tpi(0) — ip{L) = 0, that is equiv- 
alent V trap (0) = Vtrap(L) = +00. FigfJJa) shows the 
eigenspectrum for a non-optimized potential, correspond- 
ing to equally spaced potential wells (FigQJc)). The oc- 
cupied states colored in blue and unoccupied states are 
shown in magenta. Note, that the trapping potential is 
scaled down by the factor of 10 for better visibility. 

First we perform the optimization for relatively small 
level broadening 7 = 0.01-En and temperature T = OK. 
The optimized heterostructure has the eigenspectrum 
shown in FigOJb), and FigfTJd) shows the optimized 
trapping potential and the corresponding ground state 
electron density. Note, that the optimized potential has 
a dramatically decreased energy gap A£jv,jv+i between 
the highest occupied and the lowest unoccupied energy 
levels. The initial energy gap is Ai^jv+i ~ O.ll-Eo, 
and after the optimization ASjv,jv+i ~ 0.005i?n- For 
the field enhancement estimate we only need the term 
in Eq.([T]) describing the transition between levels N and 
N + 1. The decrease of the energy gap leads to an en- 
hancement of the N — > iV+1 transition by approximately 
I o° 005 — 'o °oi 2 ] ~ and the corresponding field intensity 
(and stored energy W) by a factor of w 100, compared 
to the non-optimized potential FigQJe). Potentially even 
further enhancement of the polarization may be achieved, 
if there are several unoccupied energy levels (i.e. a higher 
density of states), close to the highest occupied one. 

Since the energy W stored in the heterostructure is 
increased by approximately factor of 450 in the static 
limit (see Figfljf)), it is clear, that the vanishing of the 
gap is a significant, but not the only mechanism for the 
field enhancement. The enhancement of the transition 
dipole matrix elements, including for N — > ./V + 1 is an- 
other mechanism. However, for narrow energy levels this 
mechanism gives a relatively less significant contribution 
of a factor of w 4. 

Upon increasing the level broadening the dipole mech- 
anism becomes more significant. In Fig|2] we consider 
optimization of the trapping potential in a case of larger 
level broadening 7 — O.li^o- First, for the optimized 
heterostructure in Figl2^a) the electron density has the 
shape with two maxima, resulting in a bigger dipole re- 
sponse of the system. Second, for the optimized trapping 
potential the enhancement of the energy stored in the 
heterostructure Figj2£b) is much more modest, compared 
to Figrjjf). The optimized design provides improvement 
of the stored energy by a factor of 2.71. Third, the en- 
ergy gap AEn,n+i ~ 0.02i?o for the optimized potential 
is smaller, than in the non-optimized one (see FigQJa)), 
but it is still k 4 times larger than in FigQJb). Thus, 
the contribution due to the resonant effect and the con- 
tribution due to the enhancement of the dipole matrix 
elements are much closer, than in the previous example, 
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FIG. 2: (a) The effective trapping potential (black dashed 
line) and the corresponding ground state electron density (red 
solid line) for the optimized heterostructure. (b) Frequency 
dependence of the induced field energy in the nanostructure. 
Note, T = 0, 7 = 0.1£ - 



with narrower energy levels. Note, that the optimized 
potential for 7 = O.OlEo produces significantly less po- 
larization than the optimized potential in Figj2ja). We 
conclude, that level broadening makes the previously op- 
timized design of the heterostructure non-optimal. 

We also studied the optimization of the trapping po- 
tential at T > OK in FigJ3] Surprisingly, at the relatively 
high temperature ksT = 0.5E (here ks is the Boltz- 
mann constant), the optimization gives Vt rap (x) = as 
the best result, corresponding to a square well potential. 
A reasonable explanation of this surprising result is that 
the square well potential gives the fastest scaling of the 
energy levels En oc jV 2 , producing the maximum energy 
gap A-EjVjiV+i; and as a result, the maximum population 
difference f{E N+1 ) - f{E N ). 

In conclusion, we have studied the principles of optimal 
design of heterostructures, which function as nanoscale 
capacitors. Three main strategies to maximize the stored 
electrostatic energy were identified. The first one is the 
minimization of the energy gap AEn,n+i, the second 
is the maximization of the transition dipole matrix ele- 
ments for the electrons in the heterostructure, and the 
third is the maximization of the population difference 
/(.E/v+i) — I{En) in the system. The effectiveness of 
each strategy depends the level broadening and temper- 
ature of the system, which we treated as independent 
parameters. Note, that the first and the third strategies 
are mutually exclusive. 

Some caution is necessary, as the minimization the en- 
ergy gap between the highest occupied and lowest unoc- 
cupied levels is equivalent to making system "more metal- 
lic" . In this case the conductivity in such system should 



grow in the same way as the susceptibility as a function 
of the energy gap (AEnn+i)~ 2 , which is consistent with 
experimental findings [2j • Since for energy storage appli- 
cations one needs nanocapacitors with relatively small 
inner charge leakage, the scaling of the electrical conduc- 
tivity and tunneling effects have to be taken into account. 
Another effect, which we can not take into account within 
linear response theory is the dielectric breakdown. The 
increase of the breakdown voltage is favorable for en- 
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FIG. 3: (a) The effective trapping potential (black dashed 
line) and the corresponding ground state electron density (red 
solid line) for the optimized configuration, (b) Frequency de- 
pendence of the induced field energy in the nanostructure. 
Note, k B T = 0.5E , 7 = 0.0LE o . 



ergy storage applications, since the electrostatic stored 
energy W oc V^ ax , where the maximum work voltage 
Vmax cannot exceed the breakdown limit. In general, 
smaller energy gap will result in a smaller critical break- 
down voltage for bulk materials, however, for spatially 
inhomogeneous nanostructures more systematic ab initio 
simulations are desirable. 
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